A glass-like behavior in the low-temperature specific heat 
is a natural property of any real crystal 
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We provide a rigorous calculation of the free energy of a non-metallic crystal containing a small 
concentration of defects. The low-temperature leading contribution is found to be oc T^. This 
further gives a linear-in-T low-temperature specific heat as that exhibited by glasses. These results 
also show that, similarly to what happens in glasses, the long-wavelength spectrum of a nearly 
perfect crystal does not suffice to determine its low-temperature behavior. 



It is well recognized that the presence of defects may 
have a strong influence on the dynamics of the crystal 
lattice. The main effects come from the changes of the 
parameters characterizing the crystal and the changes in 
the corresponding structure induced by the defects. The 
former are normally modeled as local inhomogcneitics, 
e.g., local changes in the elastic moduli and/or mass den- 
sities, while the later are encoded in the form of local 
fields [l|, [3, [1, 13, H- Whatever the facet of the influence 
of defects one considers, forced oscillations of the crystal 
lattice are inevitably accompanied by the (rigid) motion 
of defects. And this motion of defects further gives rise 
to the radiation of long-wavelength acoustic waves. This 
process is behind, for instance, the dielectric losses pre- 
dicted as a result of the presence of charged defects in 
ionic crystals The aim of this Letter is to show that 
this type of processes involving defects modeled as local 
fields plays a key role in the low-temperature behavior 
of any real crystal since any real defect has a local field 
component. Specifically we show that the free energy 
of the crystal lattice has a contribution oc resulting 
from the presence of defects 0] , which further implies a 
linear-in-T specific heat at low temperatures. 

In the presence of defects, there are local couplings 
between the modes of the perfect lattice that in gen- 
eral modify the corresponding density of states. These 
changes, however, are normally quite difficult to com- 
pute. So, to reveal the influence of defects, it is more 
convenient to follow a different approach. It is under- 
stood since long ago that the infiuence of defects can 
be revealed (at least partially) by considering an aver- 
age over the defect configurations [3, H, U, i8|] . Thus the 
influence of defects is encoded, among other quantities, 
into the damping of the elementary excitations of the av- 
eraged system. Then, the potential importance for the 
low-temperature thermodynamics of the presence of de- 
fects can be clearly foreseen by realizing that they provide 
a mechanism for dissipation. In fact, through the afore- 
mentioned local couplings, the energy of lattice modes of 
small wavelengths blur into large- wavelengths ones. 

Once the elementary excitations of the system are sub- 
ject to damping, their thermodynamic properties change 
qualitatively. For example, the free energy of a damped 



oscillator has a dependence at low temperatures 



F : 



(1) 



instead of usual dependence of the undamped case . 
Here Fq is the ground-state energy, 7 is the damping co- 
efficient and ujQ is the natural frequency of the oscillator 
(here and hereafter we use units such that h = ks = !)■ 
In the view of Eq. ([T]) it is clear that there will be, e.g., a 
linear-in-T contribution to the specific heat due to these 
damped oscillations, as already was argued in [ll|, [T^ . 
Similar speculations were made in the realm of amor- 
phous solids These papers deal with the damping 
phenomenologically and, in the case of amorphous solids, 
a peculiar dependence on temperature was just assumed. 
In the present Letter we address nearly perfect crystals 
and, making no assumptions about the damping, we de- 
rive these results from rigorous microscopic calculations 
showing that they are related to the motion of defects. 

To understand that there exists an unavoidable mo- 
tion of defects that, even at zero temperature, causes the 
damping of phonons, let us consider an acoustic standing 
wave in the presence of a defect. For the sake of concrete- 
ness, we will restrict ourselves to point defects acting as 
dilatation centers 0] ■ The standing wave should be imag- 
ined as a forced oscillation of the crystal whose frequency 
Lu and wavevector k are consequently independent one 
each other. This allows us to compute the zero-frequency 
limit of the damping constant for this oscillation, what 
is precisely the constant 7 that enters Eq. ([T|). As re- 
gards the defect, the gradients in pressure associated with 
the standing wave induce forces that make it moves [iJl ■ 
These forces vary in time with the frequency u and, once 
the defect is in motion, this represents a source of acous- 
tic radiation. The frequency of this radiation is again w, 
while its characteristic wavevector is w/c (c is the cor- 
responding velocity of sound). Notice that, as we are 
considering the limit w — > 0, this later wavevector lu/c 
is much more smaller than the wavevector k of standing 
wave. Already known formulas of acoustics tell us that 
the intensity of the acoustic radiation is oc p^ . This 
energy evidently comes from the wave inducing the mo- 
tion of the defect. Putting it differently, this wave losses 
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its energy and therefore is subject to damping. And this 
damping is just a frequency-independent (Ohmic) damp- 
ing as in Eq. ([T]) It is worth mentioning that the 
acoustic degrees of freedom of small wavevectors then 
play the role of a reservoir that accepts this energy, what 
is in tune with the system-subsystem picture implicit in 
Eq. (P) (see Ref. 0). 

Let us proceed with more detailed calculations. To 
take into account the presence of dilatation centers we 
introduce the term 



-huii 



(2) 



in the Hamiltonian density of the system, where Uu is the 
dilatation of the system, i.e., the relative volume change, 
and 



(3) 



with ha being a constant characterizing the change in 
the crystal volume caused by the n-th defect. As we 
have mentioned, the key point in our considerations is the 
accounting for the fact that the positions of the defects 
change as 



u 



(4) 



due to the oscillations of the lattice, where u is the dis- 
placement vector describing these oscillations. 

We are interested in computing the temperature de- 
pendence of the free energy resulting from the thermal 
motion of the lattice in the presence of defects. Mat- 
subara's method provides us a convenient diagrammatic 
scheme, within which the thermodynamic properties are 
obtained from so-called temperature Green's functions 



Here denotes the ordering with respect to the (imag- 
inary time) variable r e [0, 1/T], Ui and Ui are Matsub- 
ara operators in the "interaction representation" defined 
from the Schrodinger displacement operators iii and m| 
as 



F.(r,T) = ^e^"-<(r)e-^°^ 



(6) 
(7) 



where p is the density, and ^ = ^(1/T)is the Matsub- 
ara c5^-matrix defined as 



^(r) = Tr exp 



i?dcf(r')rfr' 



(8) 



where Hdcf{T) = e^^'^H^ci e ^"^ ■ In deriving these for- 
mulas, the total Hamiltonian is separated into Ho + H^cf, 
where Hq is the Hamiltonian of the ideal system and Hdd 



is due to the defects (we assume that the system is har- 
monic, so the interaction is solely due to the defects). 
Averaging is carried out with respect the states of the 
ideal system. 

The correction to free energy associated with the pres- 
ence of defects can be obtained from connected diagrams 
contributing to the ^-matrix as AF = — r((=5^)conn — 1) 
Q. Thus, the basic ingredients in this perturbativc ap- 
proach are the zero-order temperature Green's functions 
and the defect vertices. 

In Fourier space, the temperature Green's functions for 
the acoustic phonons are 



(9) 



in the ideal isotropic case [l7|. Here a;„ = 27rnT are 
Matsubara frequencies [n = 0, ±1, ±2, . . . ) and LOi{k) = 
cik and LjJt{k) = Ctk are the energies of the longitudinal 
and transverse phonons, respectively. 

To find out the factors associated with the defect ver- 
tices, let us consider the contribution to the defect Hamil- 
tonian due to a given defect. In accordance with Eq. ([2]), 
this can be written as 



dct.O 



[k . u(k)], 



^-ik-(ro-(-u) 



(10) 



where u(k) is the Fourier transform of the acoustic dis- 
placement vector: 



u(r) 



dk 

(2^ 



.u(k)e 



— ik] 



(11) 



This vector can be split into u = u° -I- u', where u' rep- 
resents the displacement with respect to the (static) po- 
sition of equilibrium of the defect, Rq = Tq + u°, due to 
the oscillations of the lattice. Furthermore, by expand- 
ing Hdci.o in terms of u', we keep the nontrivial terms of 
lowest order: 



ho J ^[k.u'(k)][k.u'(Ro)]e 
J (^[k-u°(k)][k.u'(R 



-?k-R,o 



(12) 
(13) 



In the second term we have the Fourier transform of the 
equilibrium dilatation of the crystal — ik • u°(k). To the 
lowest order in ha, this quantity is such that 

ik- u°(k) = hk,G^j{k)Y,hae'^"^\ (14) 



where G^j{k) = p-%^p{k,0). Furthe rmore, taking into 
account Eq. PT]) . Eqs. and can be rewritten as 

I f dk f dk' 



2 J (27r)3 J (27r)3 
I f dk f dk' 



2 J (27r)3 J (27r)3 



h,,.^oiKk'K{k)u'Jk'), (15) 



5.,,o(k,k')i4(k)u;(k'). (16) 
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Here 



with 



/ii,-o(k,k') = hoihk, + fc^fc')e-*(''+'^')-^«, (17) 



5ij,o(k,k') = ^^jj oae" 



i(k:+k')-Ro 



(18) 



where 



g^J,Oa = -hoha J ^^^fc, fcj fc, fc^G,™ (k)e 

(19) 

The terms (fTS)) and p6)) give rise to two types of ver- 
tices in the diagrams contributing to the free energy, as 
shown in Fig. [TJ In these diagrams each hne corresponds 
to (h, ujn) while defect vertices indicated by crosses 
are associated with the factor 

p-i/iy-,(k,k')5..„.^„,, (20) 
and defect vertices indicated by circles are associated 



'ig,,,a(k,k')5. 



(21) 



where (k, Wn) and (k',u;„') are the arguments of the 
Green's functions meeting the corresponding vertex. At 
these vertices, sum with respect to the subindices of the 
Green's functions is also carried out taking into account 
the subindices in the above factors. 

For the sake of simplicity, let us assume that defects 
do not change the total volume of the system and that 
there is no correlation between them. In this case the 
coefficients ha can be taken as ±/i, with the same con- 
centration of defects for each sign N± — N/2. 

The contribution of diagrams containing an odd num- 
ber of cross vertices then vanish after averaging over the 
defect configuration. So the lowest order correction given 
by these diagrams is due to those of the type (2a) in Fig. 
[T] Furthermore, these diagrams give a nonzero contribu- 
tion after averaging only if the crosses refer to the same 
defect. As a result we find a contribution proportional to 



J 



Nh'^T 



dk 



dk' 



(27r)3 J (27r) 



(^kikjkkki + 2kikjki^ki 



^°^(k',c.„), 



(22) 



r 



where h' ~ h/p. At low temperatures, the most im- 
portant term comes from the first product of fc's in the 
integrand. This term reduces to 



(2^)3 i (2^)3 



(2^)3 7 (2^)3c.2+^f(fc) 



(23) 



Here we can write 
T 



E 



uji{k) - UJi{k') 



2uji{k)u;i{k')[ujf{k)^ujf{k')] 



n[uji{k)] n[LUi{k')] 



Ljfik')~u;fik) V Wi(fc) ujiik') 

[,.,2 _i_ ,,,2rz,Mr,,,2 _i_ ,.,2fl,l\]l,.,2 j_,.,2fl,!\ 



[o:l + u;fik)][u;l+cof{k')][u;l+u;Hk')r 

(24) 

where n(aj) = [cxp{lli/T) — 1]^^ is the Bose-Einstein dis- 
tribution function. The first term at the rhs of p4l) docs 
not depend on the temperature. At low temperatures, 
the contribution of the second term can be estimated as 



-2Nh'- 



dk r dk' k^kln[uji{k')] 
(2^7 (27r)3 u;f{k)cji{k') 



(25) 



where fcmax is the radius of the Brillouin zone. One can 
see that both the remaining sum in Eq. and the term 
(fTH)) give similar contributions ^ to the free energy. In 
the later case, a nonzero contribution is already obtained 
from diagrams of type (2b) in Fig. [T] [notice that the 
factor (PT|) already is ~ h"^]. 

At low enough temperatures, the contribution we 
have obtained due to the presence of defects can surpass 
the one of the ideal lattice (oc T'^, see e.g. Ref. jlG|'). 
Eq. (|25p indeed coincides with the free energy estimated 
for a set of harmonic oscillators, with the dispersion law 
cj;(fc) = cik and a damping constant 7 ^ Nh''^kf^^^/ci 
[see Eq. ([TJ]. As wc have argued, this coincidence is not 
accidental. The phonon lifetime indeed can be computed 
within a pcrturbativc approach similar to employed here. 
As a result one finds that, for large enough wavevec- 
tors, acoustic phonons are precisely subject to an Ohmic 
damping due to the motion of defects [ifll- The reason- 
ings given in Refs. [H!, are therefore apphcable when 
estimating the crossover temperature. This temperature 
is found to be attainable experimentally (up to 1 K for 
moderate concentrations of defects). 

It is worth mentioning that, as a result of higher order 
terms in the above expansion of the c5^-matrix, one sim- 
ply obtains corrections to the coefficient of the term 
in the free energy (also new terms with higher powers 
of T, unimportant at low enough temperatures). When 
it comes to the comparison between theory and experi- 
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(1a) (2a) (2b) (3a) (3b) (4a) (4b) (4c) 
FIG. 1: Diagrams contributing to the free energy. 



ments, this implies that the lowest-order result we have 
found is robust (qualitative agreement is expected). But 
it also have a more profound reading. The fact that the 
density of defects can be increased without altering the 
asymptotic low-temperature behavior allows one to spec- 
ulate that the result is valid even in a glass limit (i.e., for 
large concentrations of defects). Indeed, already with a 
small concentration of defects we obtain, e.g., a linear- 
in-T low-temperature specific heat as that exhibited by 
glasses (see e.g. Ref. [l^). This linear- in-T specific heat 
is normally considered as a fingerprint of glasses. How- 
ever, we see that these distinctive features can already be 
found within the conventional scenario of a crystal with 
a not necessarily large concentration of defects. 

We have so far restricted ourselves to reveal the influ- 
ence of defects through the local dilatation they induce. 
Real defects, however, have a broader influence Any 
impurity rarely acts as a dilatation center exclusively. On 
the contrary, they generally induce complicated force dis- 
tributions. Besides defects may be individually charged 
and/or they can act as local electric multipoles and, con- 
sequenly, they will move due to electric fields varying 
in time (even inhomogeneous ones). Beyond that, the 
properties of the corresponding crystal can be affected 
by defects as more general multipoles, say "optic mul- 
tipoles" 0, [3, that also will oscillate if there exists a 
wave of the corresponding optic coordinate. So the ele- 



mentary excitations of the system inevitably will induce 
the motion of the defects. And this motion will gener- 
ate a long-wavelcngth acoustic radiation. The process on 
which our results base is therefore universal in the sense 
that takes place for any type of defect. Consequently, the 
same thermodynamics is expected for any real crystal. 

Let us mention, in closing, that the contribution we 
have found associated with the motion of defects is not 
a mere correction of previous results. It indeed may be 
the leading contribution by virtue of the low tempera- 
ture of the system. Moreover, contrary to what happens 
in ideal case, lar ge w avevectors are involved in this lead- 
ing contribution [19| . This represents a deep change with 
respect to the paradigm for the low-temperature behav- 
iors (see e.g. Ref. H, [l^l), according to which these 
behaviors are associated with the macroscopic proper- 
ties of the corresponding system. In some sense, it takes 
place the same situation as in glasses. That is, universal 
low-temperature behaviors are found but cannot be un- 
derstood as uniquely due to the macroscopic properties 
of the system. 
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